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Role of unstable periodic orbits in phase transitions of coupled map lattices 
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The thermodynamic formahsm for dynamical systems with many degrees of freedom is extended 
to deal with time averages and fluctuations of some macroscopic quantity along typical orbits, and 
applied to coupled map lattices exhibiting phase transitions. Thereby, it turns out that a seed of 
phase transition is embedded as an anomalous distribution of unstable periodic orbits, which appears 
as a so-called q-phase transition in the spatio-temporal configuration space. This intimate relation 
between phase transitions and q-phase transitions leads to one natural way of defining transitions 
and their order in extended chaotic systems. Furthermore, a basis is obtained on which we can treat 
locally introduced control parameters as macroscopic "temperature" in some cases involved with 
phase transitions. 
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I. INTRODUCTION 

For the past few decades chaotic dynamical systems 
with a few degrees of freedom (DOFs) have been investi- 
gated theoretically, numerically, and experimentally with 
enthusiasm, which has brought various insights about 
them. Since one cannot follow individual trajectories 
in chaotic systems by any means, one of the subjects 
attracting interest is evaluation of dynamical averages, 
namely asymptotic time averages and fluctuations of 
some observables along typical orbits. The thermody- 
namic formalisnillL 4] , which is frequently used for multi- 
fractal analysis E 1 1 is exploited for this purpose [4] and 
the concept of dynamical averaging has been remarkably 
developed by means of unstable periodic orbit expan- 
sion, trace formulae, and dynamical zeta function, which 
reveal the role of unstable periodic orbits (UPOs) as a 
skeleton of chaos On the other hand, the thermody- 
namic formalism is sometimes discussed in the context 
of phase transitions, called q-phase transitions. This is 
not a transition dealt with in statistical mechanics which 
involves large fluctuations of thermodynamic quantities 
and occurs only in the thermodynamic limit, but a tran- 
sition with large dynamical fluctuations of observables 
which occurs in the long-time limit. It has been shown 
that the large fluctuations reflect the dynamics and q- 
phase transitions indicate a singular local structure of 
the chaotic attractor, such as homoclinic tangencies of 
stable and unstable manifolds and band crises [6| . 

Despite the understanding of low-dimensional chaotic 
systems, less is known about spatially extended systems 
whose number of active DOFs is large or inflnite. This 
is partly because of difficulty in treating concepts such 
as measures for infinite-dimensional dynamical systems 
in a mathematically proper way 0, and partly because 
of computational complexity; e.g. with regarding to the 
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UPO expansion, not only does the number of UPOs grow 
exponentially with increasing DOFs, but even finding one 
UPO becomes much more laboring. However, the num- 
ber of DOFs one can numerically investigate increases 
gradually, which makes it possible for various theoreti- 
cal concepts and methods to be extended and applied to 
high-dimensional chaos [1, Q . It leads to reveal several 
suggestive properties intrinsic to spatially extended sys- 
tems, which have been reported recently. For example, it 
was found that one can reproduce macroscopic quantities 
of turbulence only from a single UPO (Tol. Ill|. 

One of the most striking manifestations of high dimen- 
sionality is the occurrence of phase transitions. In the 
case of coupled map lattices (CMLs), i.e. lattices of inter- 
acting dynamical systems whose time evolution is defined 
by a map, logistic CMLs are known to display nontriv- 
ial collective behavior which cannot be observed in equi- 
librium systems, and transitions between two types of 
collective behavior can be regarded as phase transitions 
[1 , [1^ . Another interesting example of non-equilibrium 
phase transitions is 2-dimensional CMLs which exhibit a 
continuous phase transition similar to that of the Ising 
model [3, TTsj . The existence of a new universality 
class was numerically shown for such CMLs with syn- 
chronously updating rules, while those asynchronously 
updated belong to the 2D Ising universality class [15| . 
Recent studies suggest that the Ising-like transitions of 
synchronous CMLs and the onset of the non-trivial col- 
lective behavior of a logistic CML bclongto the same 
universality class, i.e. the non-Ising class [13|. 

Although many interesting properties of non- 
equilibrium phase transitions have been found out, 
there seems nevertheless no consensus on the usage 
of the term "phase transition" in dynamical systems. 
Theoretically, it can be defined as a qualitative change 
in the statistical behavior of typical orbits in a single 
mixing attractor which does not change topologically 
[il,!j9], by which we exclude bifurcations coming 
up even in finite-dimensional dynamical systems. For the 
definition of "qualitative change," the analogy with that 
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in equilibrium phase transitions is used. There are two 
complementary manners of characterizing equilibrium 
transitions [20| : one is after Ehrenfest, where n-th 
order phase transitions are identified as divergence 
or discontinuity of some n-th derivative of the free 
energy. The other is after Gibbs, where first-order phase 
transitions correspond to a change in the number of the 
pure Gibbs measures, or macrostates. Analogues of the 
latter have been adopted in the context of dynamical 
systems since no free energy appears useful: if we 
consider interaction in a formal Hamiltonian on the 
space-time configuration space and define a free energy 
from it, then the analyticity of the free energy is a 
very delicate problem [3, [2l| and too complicated to 
relate to phase transitions. On the other hand, if we 
define a free energy from purely probabilistic measure 
approach as we will see below, then it is identically zero 
and thus analytic in the whole parameter region, even 
at criticality, due to a strong constraint which comes 
from a normalization of the measure [l^. Therefore, 
transitions have been defined in the Gibbsian sense 
, Il9l | , that is via a change of the number of natural 
measures, which corresponds to first order transitions. 
This definition, however, cannot characterize higher 
order transitions as definitively, so we have to make 
use of more subtle phenomena, such as spontaneous 
symmetry breaking, divergence of a correlation length, 
formation of an infinite cluster, and so on. It is true that 
they are closely related to phase transitions, but would 
not prescribe them as quantitatively as equilibrium 
counterparts do. Thus, it is desirable to develop another 
way to characterize phase transitions in extended chaos, 
including higher order ones. 

Another issue involved with phase transitions in ex- 
tended chaotic systems is absence of macroscopic "tem- 
perature" which controls the systems. Some locally de- 
fined parameters such as coupling strength and diffusion 
coefficient have been used as ad hoc substitutes for tem- 
perature (e.g. [13, [l^), while its theoretical grounds 
remain to be clarified. This treatment is based on an 
assumption that such local parameters are direct barom- 
eters of macroscopic properties. This is however not at all 
trivial, as we can see for example from studies of effective 
temperature in non-equilibrium systems [2^. Although 
we can argue the issue to some extent by renormalization 
group approach, it cannot deal with concrete systems. 
Therefore, it is desirable to have a basis on which we can 
connect locally defined time evolution rule of a system to 
macroscopic properties. 

In the present paper, we characterize phase transitions 
in extended chaotic systems, namely CMLs, including 
both equilibrium and non-equilibrium ones. The peri- 
odic orbit expansion and the thermodynamic formalism 
are adapted for such systems, by which the relation be- 
tween q-phase transitions and (actual) phase transitions 
is investigated. The main outcomes are that (1) a rather 
quantitative way to define non-equilibrium phase transi- 
tions with their order in the Ehrenfest's sense is proposed. 



and (2) a basis is obtained on which microscopic control 
parameters can be handled similarly to temperature in 
some cases involved with phase transitions. Note that 
it is not the aim of this paper to give a mathematically 
rigorous argument, which is often highly delicate in this 
field 0, ISH and may limit an attainable conclusion. In- 
stead, we shall devote ourselves to obtain a physically 
plausible picture. 

This paper is organized as follows : we first review 
the idea of UPO ensemble (Sec. UlAl . and on its 

basis the thermodynamic formalism is formulated to deal 
with dynamical averages and fluctuations of some macro- 
scopic quantity in chaotic systems with many DOFs (Sec. 
IIIBI) . A corresponding partition function and topological 
pressure, or "free energy" , are defined and the moments 
are obtained as its differential coefficients. Then we ap- 
ply it to a ID Bernoulh CML (Sec. which can be 
regarded as a deterministic model of the ID Ising model 
as is summarized in Sec. IIII Al After we mention the 
computational procedure for the thermodynamic formal- 
ism (Sec. IIII Bp . we show that an anomalous distribution 
of UPOs exists in such a system with phase transitions 
(marginal transitions in this example), which can be re- 
garded as a seed of the Ising transition (Sec. IIII C|) . The 
seed is embodied as a q-phase transition. Another exam- 
ple is a ID repelling CML which exhibits a non-marginal 
transition (Sec. IIV[) . This is a solvable case, hence we can 
explicitly see the relation between phase transitions and 
q-phase transitions. Section |V] is assigned to the discus- 
sion and conclusion. Note that the terminologies "phase 
transition" and "q-phase transition" are specifically dis- 
criminated throughout this paper. 



II. THERMODYNAMIC FORMALISM 
FOR EXTENDED SYSTEMS 

A. UPO ensemble 

First, we review the concept of UPO ensemble [lll[2^. 
on which the following thermodynamic formalism is 
based. This and the next subsections are assigned to 
show the grounds for our arguments in the rest of sec- 
tions and the range of their applications. 

Consider a dynamical system with discrete time, 
^t+i _ a;t = [xq,x{,--- ,a;^_i], where N de- 

notes the number of DOFs and is large. Our goal for 
the time being is to obtain the dynamical average of 
an arbitrary macroscopic observable A{x), which is de- 
fined as a function of the dynamical variable x. Here 
the term "macroscopic observable" represents a quantity 
obtained by taking the average over the DOFs of the 
system. Suppose the system is ergodic, the long-time av- 
erage (A)time = lim„_^oo(l/n) J2t=o ^(a;*) is equal to the 
phase space average (A)^ = / A{x)^{dx) for almost all 
initial conditions x^, where fj, denotes the SRB measure, 
or the natural invariant measure. For mixing and hyper- 
bolic systems, the following relation between the natural 
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invariant measure of a subset S and UPOs holds \24 



B. Thermodynamic formalism 



fi{S) ^ lim V 



(1) 



{[x]-xOes} 



Here, [x] = x'^x^ ■ ■ ■ x^^^ indicates an UPO of period p 
and therefore the sum in Eq. ([1]) is taken over aU the 
UPOs of period p which start from S and return to it. 
A([a;]) is a positive Lyapunov exponent per 1 DOF 



1 

1=0 



(2) 



where {A^([a;])} denotes a set of positive exponents of 
the UPO [a;]. Note that we sometimes cah A([a;]) simply 
"Lyapunov exponent" as long as it does not cause any 
confusion. From Eq. ([1]), q^p^^H^]) can be regarded as 
the probability measure of the UPO ensemble. 

For fixed p longer than the time required for mix- 
ing, denoted by Tmixing, Eq. ([T]) holds approximately, 
so that {A)fj_ can be estimated from the ensem- 
ble average {A{x"))upo = ^2^^^ A{x°)e-P'^^<^^'^^^ ~ 
J A{x'^)ii{dx'^) — {A)fj_. Moreover, since A([a;]) is in- 
variant under cyclic permutation of [a;] = x^x^ ■ ■ ■ x^^^, 
{A{x*))upQ ~ {A)f^ also stands, and consequently for the 
average along an UPO 



p-i 



Ai[x])^-Y,A{x 



t=o 



the following relation holds if p ^ ^mixing 

{A{[x]))^PO ^ 5] A([a.])e-^^^(W) ^ (^a), 

[CD] 



(3) 



(4) 



Note that the period p required to make Eq. ^ con- 
verge can be shorter than Tmixing if N is sufficiently large, 
thanks to the law of large numbers [1^ . The estimation 
of {A)fi from A([a;]) defined by Eq. ([3]) is preferable to 
that from A{x*), because its variance 



a(A(M))^PO^((^(W)-(A(M))upo) 



UPO 



p-1 



UPO 



(5) 



r=0 



does not exceed the variance of A{x*), a{A{x*))'^pQ = 
{A{x*) , A{x*)) ~ o'(^)^- If we assume that 
the autocorrelation function decays exponentially 
(yl(a;*), A(a;*+^))uPQ e"^/^" with the correlation 
time 1 ^ To <C p/2, the ratio of the 2 variances is 
aiAi[x]))lp^ {2ro/p)a{A)l. 

As is shown above, the phase space average of the 
macroscopic quantity A(x) and the lower bound of its 
fluctuation are obtained from the UPO ensemble treat- 
ment, which are approximately equal to those time aver- 
aged along typical orbits. 



In this subsection, we introduce an appropriate parti- 
tion function to deal with dynamical averages and fluc- 
tuations in extended systems, that is 



-pN[qX([x])+l3A([x])] 



(6) 



where the sum is taken over all of the UPOs whose pe- 
riod is p. The summation without the second term in 
the exponential represents the Lyapunov partition func- 
tion ^26]. Variables q and P inserted in Eq. ^ are aux- 
iliary ones, which can be regarded as inverse tempera- 
ture mathematically, but of no particular physical sig- 
nificance. However, since we can change the dominant 
terms in the sum of Eq. © by varying q and f3, they play 
essential roles in the following argument. The real sys- 
tem corresponds to {q,f3) = (1,0), where the summands 
in the partition function coincide with the probability 
measures of the UPO ensemble, hence we call it physical 
situation hereafter. Note that the partition function ^ 
is similar to that introduced by Fujisaka and Inoue [J], 
but here we explicitly consider the scaling dependence on 
the number of DOFs N as well as the period p in order 
to argue phase transitions. 

The relation to the space-time Gibbs measure should 
be also referred to. The space-time measure is often in- 
troduced as a measure of refinement elements (so-called 
cylinder) under the symbolization fl^, [l3| and the ac- 
companying partition function is equal to that in Eq. ([B]) 
with (g,/3) = (1,0). The partition function mentioned in 
this paper is constituted by adding an observable A([a;]) 
to the argument of the exponential term and introducing 
"temperature" parameters based on the thermodynamic 
formalism. What is essential in the concept of the space- 
time measure is that we can consider its configuration 
space to be the {d -\- l)-dimensional space-time compris- 
ing the d-dimensional space and the 1-dimensional time, 
which remains valid after the extension. That is to say, 
the partition function defined here are natural extensions 
of those of the space-time measure, which means we can 
exploit plentiful knowledge in the equilibrium statistical 
mechanics to the problem of spatio-temporal chaos. 

The corresponding "free energy" , which is called the 
topological pressure in the context of dynamical systems 
P, Q but we call it here the generalized Massieu function 
(GMF), is defined by 



(7) 



Note that the sign of Eq. ([7]) is set opposite to the con- 
ventional definition of the topological pressure so as to 
mention a minimum principle of it. Since both A([a;]) 
and A([a;]) can be regarded as intensive densities per 
unit time and one DOF, they remain finite in the limit 
p,N ^ oo and the GMF ^'(g, (3) is expected to converge 
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in that limit. Especially, since the measure of the whole 
phase space is one, Eqs. ([1]), ©, and ([7|) yield 



lim ^'(1,0) = 0. 



(8) 



This constraint must be satisfied at the physical situa- 
tion regardless of values of control parameters. This fact 
prevents us from defining phase transitions just by the 
sigularity of the "free energy" with respect to parame- 
ters. We shall see, however, that by introducing a gener- 
alized probability measure in Eq. we make room for 
the singularity with respect to q and (3, which is called q- 
phase transition ^] , and thus we are in fact able to relate 
actual phase transitions (with respect to parameters) to 
the singularity of the "free energy." This point will be 
clarified in Sec. IIII CI 

The ensemble average and fluctuation of ^([a;]) defined 
by Eqs. @ and (O respectively are obtained from the 
differential coefficients of the GMF by 



(A),,^(A([a3]))uPo = ||(l,0), 



a{A)l > aiA{[x]))l^o 



1 6)2^ 



(9) 



(1,0), (10) 



where the UPO average (A([a;]))upo is redefined as 
(^(W))uPO ^ Ew^(W)e-^^'"=^lVEwe-^^'"^» in 
order to moderate errors due to the finite-size effect. 
These relations are completely analogous to counterparts 
of the canonical statistical mechanics and therefore all 
moments of ^([a;]) can be obtained by differentiating the 
GMF up to the requisite order. The average and variance 
of the positive finite-time Lyapunov exponent per 1 DOF 
can be also acquired without replacing the definition of 
A{x) by them; 

(A([a.]))^c.(A([a;]))uPo = 1^(1,0), (11) 

a{\{[x]))l a{\{[x]))l^o = -^^(1' «)■ (12) 

The positive (infinite-time) Lyapunov exponent per 1 
DOF can be obtained by taking a limit p ^ oo in 
Eq. pT|) . Moreover, the equalities and inequality ([9])- 
(fT2|) hold precisely in that limit. They are expected to 
be good approximations for a finite period p, at least if 
P ^ ^mixing, as is mentioned in the previous subsection. 

Now we consider the statistics of UPOs, namely the 
distribution of UPOs with respect to their macroscopic 
properties. Let ri(A, A)dAdA denote the number of UPOs 
whose positive Lyapunov exponent A([a;]) and macro- 
scopic quantity ^([a;]) are within the range of A-A + dA 
and A-A -I- d^ respectively. Suppose the system is ho- 
mogeneous, in other words the system consists of iden- 
tical DOFs and thus it can be viewed as an ensemble of 
smaller coupled subsystems, we can assume the following 
functional form of A)d\<lA: 

n{\ A)d\dA - p{\, A; p, iV)eP^^(^^^MAdA. (13) 



Here p{\,A]p,N) is a "coefficient" into which all factors 
are pushed whose dependence on pN is not exponential. 
H{\, A) is a concave function, which is considered to be 
a topological entropy per 1 DOF under the restriction of 
A([a;]) e [A, A + dA] and A{[x\) e [A,A-hdA]. Roughly 
speaking, the expression is justified by the large de- 
viation theorem because both A ([a;]) and ^([a;]) can be 
regarded as the averages over pN variables which corre- 
late to each other with a specific correlation time and 
length. By making use of the distribution function ([Te 
to calculate the partition function ([H]), we obtain 



n{\A)dMA 



p(A, A- p, iV)e-P^[«^+^^-^(^^^)ldAdA. (14) 

If the product of the period of the UPOs and the number 
of the DOFs, pN, is sufficiently large, the saddle point 
approximation is applicable, that is, only the vicinity of 
the point (A, A) = {\{q,f3),A{q,(3)) where the integrand 
has a maximum contributes to the integral (jl4[) . The 
conditions imposed on X{q, P), A{q, (3) are 

dH dH ^ 

9' = 



dA 



9A2 



< 0, 



det ^ 
\dxaA 



(15a) 



> 0, (15b) 



where all differential coefficients of H{\, A) are evaluated 
at (A, A) = [\{q, (3), A{q, P)). Using the saddle point 
approximation to Eq. and substituting it for Eq. ([7]), 
we obtain 



*(<Z,/3) 



min \q\ + f3A- H(X, A)] , 

\,A 



(16) 



or 



^{q,P) ~ qX{q,/3)+/3A{q,/3)^H{X{q,/3),A{q,/3)). (17) 

These equations hold rigorously in the limit p, iV — > oo. 
Equation (|16p can be regarded as a principle of minimum 
free energy in the sense that A and A which are dominant 
in the partition function (|14p are selected out to minimize 
the corresponding GMF. The relation p?|) accompanied 
by Eq. (|15ap is the Legendre transformation and thus 



9* 

A(g,/3) -—(<?,/?), 



Aiq,f3)^^{q,P), (18) 

which are obtained by differentiating Eq. (|17|) by q or 
/3. The comparison of Eq. HHD with Eq. © and (HID 
yields the relations (A{[x]))\jpo — A{1,0), (A([a;]))uPO — 
A (1,0), which appear to be natural because the right 
hand sides represent the dominant A and A at the phys- 
ical situation. In addition, the concavity of H{X, A) 
and the relations (|15ap . (fTT]) . and ^ yield a vision of 
a general form of the function H{X,A). It is expected 
to be tangent to a plane if = A at the physical point 
{q, f3) = (1, 0), as illustrated schematically in Fig. [U and 
the tangent point represents a state which is observed 
physically. 
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tangent point (;^(1,0),A(1,0)) 



FIG. 1: A schematic view of an expected form of the topo- 
logical entropy spectrum H{X,A). 
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FIG. 2: The local Bernoulli map at a site i. The time evolu- 
tion of A- is defined in Eq. pO]! . 



III. ANALYSIS OF THE ID BERNOULLI CML 
A. Model 

The map we first analyze is a Bernoulli CML, whose 
2D version was originally proposed by Sakaguchi [23] and 
its ID version was introduced later by Kawasaki and Sasa 
[ll|. In the present work, we investigate the ID model, 
which we describe below. 

Consider a 1-dimensional lattice which consists of N 
lattice points i — 0, 1,-- - , N — 1. Dynamical variables 
{xi, Ai) E [—1, 1] X [—1, 1] are assigned to each site i, and 
in addition, a "spin" variable Si is defined as 



-1 if -1 < .Tj < A, 
-1 if A, < X, < 1. 



(19) 



. t+i 



tanh 



2(x[+s|) 



for odd 



A*+i 



A*, 



for even (i — t), 
(20) 

with periodic boundary condition = Sq. The updat- 
ing is done alternately with respect to the parity of the 
site number i, that is to say, sites with odd i are updated 
at even time t while those with even i are renewed at 
odd t. The total number of the sites N is supposed to be 
even in order that the alternately updating rule is com- 
patible with the periodic boundary condition, /(cc*; A*) 
is a Bernoulli map, illustrated in Fig. O As can be seen 
from Eq. ((20)) and Fig. [2l A* is a discrete variable and 
behaves as a dynamical parameter which describes the in- 
teraction between nearest-neighbor sites. Therefore, we 
consider only x* to be a dynamical variable and apply the 
formalism stated in the previous section. The magnitude 
and the tendency of the interaction are determined by 
the absolute value of k and its sign, respectively. For 
positive (negative) k, A* moves in the direction so as to 
make it more probable that the spin s* becomes parallel 
(antiparallel) to the neighboring spins, hence the inter- 
action is ferromagnetic (antiferromagnetic). 

The Bernoulli CML has several remarkable proper- 
ties, as demonstrated by preceding studies [HI, 113' which 
should be pointed out here. First, the dynamics can be 
expressed in terms of the symbolic dynamics with sym- 
bols s = {si}^o^. In other words, the partition {Us} of 
the phase space, whose element Us corresponds to a spin 
configuration s, is generating and thus every orbit is spec- 
ified by an infinite sequence of symbols s^s^s^ ■ ■ ■ . Espe- 
cially, note that every UFO has a one-to-one correspon- 
dence to a finite length permutation [s] = s^s^ ■ ■ ■ s'p^^. 
The most significant feature of the Bernoulli CML is that 
it respects a detailed balance and the resulting probabil- 
ity measure of a subset Us coincides with the canonical 
distribution of the ID Ising model [O, namely 



n{Us) oc exp I - ^ s.iSi+1 



(21) 



With this spin, the time evolution of (a;*. A*) is written 



Therefore the Bernoulli CML can be regarded as a deter- 
ministic model of the ID Ising (anti-)ferromagnetism in 
its equilibrium state and the interaction parameter k cor- 
responds to the inverse temperature. Since the marginal 
phase transition occurs in the ID Ising model at the zero 
temperature limit, this Bernoulli CML shows a transition 
in the strong interaction limit |fc| — s- oo. 
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B. Application of the thermodynamic formalism 

The thermodynamic formahsm in Sec. Ill Bl is made use 
of to analyze it. We adopt the Ising interaction energy 
per 1 spin for a macroscopic quantity 

1 

^(^) = -nT. ^'^'+1- (22) 



Substituting it and the Lyapunov exponent given from 
the slope of the function /(x*; A-) into Eq. we can 
obtain the following expression of the partition function 



Z. 



q,l3 



{sj.fc} 



^ ^ n.n. 



kSf,k' - I (lncoshfc)y^3j_i,fcSj,fc- 



pNq 



ln(4coshA:) 



(23) 



where a space-time configuration [s] oi N x p symbols is 
reduced to a 2-dimensional array {sj,fc} oipN/2 spins by 
exploiting the constraint s*^^ = s* for even (i — t), which 
is outlined in Fig. [31 and ^ indicates a summation 
over all pairs of neighboring spins after the spin reduc- 
tion. Equation ((23)l shows that, for positive {f3 + kq/2) 
and q, the interaction between spins comprises helical 
ferromagnetic part and spatial antiferromagnetic part. 

To calculate numerically the accompanying GMF de- 
fined by Eq. ([7]) in the limit p ^ oo with fixed N, or 
N oo with fixed p, it is well known that the zeta func- 
tion method is a powerful tool to accomplish it 0, . In 
the present analysis, however, we keep both p and N fi- 
nite in order to maintain the formal equivalency between 
space and time in Eq. P5)) and to exploit knowledge on 



t\i 123 




reduction 



(a) 



t\i 


1 23-- 







1 




2 




3 






(b) 



FIG. 3: (Color online) A schematic illustration of the spin re- 
duction from an UFO [s] to the corresponding 2-dimensional 
array {sj,fc}. (a) The spin configuration of an UPO [s]. Spins 
form clusters of length 2 in the time direction under the up- 
dating rule of Eq. (|20p . The red solid line and the blue broken 
line indicate spin interaction which comes from the Lyapunov 
exponent A([s]) and the observable A([s]), respectively, in 
Eq. (0). (b) The 2-dimensional spin array {sj^t} obtained 
by the reduction, by which each cluster is reduced to a single 
spin located at the odd (i — t). The 2 kinds of spin inter- 
action turn into helical ferromagnetic part (purple solid line) 
and spatial antiferromagnetic part (purple dotted line) sup- 
pose both (/3 -I- kq/2) and q are positive. 



the equilibrium statistical mechanics. Since the GMF 
has the same form as the Helmholtz free energy and we 
know that the GMF at physical situation is zero in sys- 
tems without escape (see Eq. ([SJ), we adopt the com- 
putational method to calculate the difference of the free 
energy [28]. From Eqs. © and (O, we obtain 



-pN[AqX([s])+Al3Ai[s])] 



(24) 



where A* = ^(g -I- Ag,/3 -f A/3) - '^{q,(3) and (• • •)g,/3 
denotes the UPO ensemble average with the probability 
distribution f{[s]) = Q-vN[q\([s])+i3A([s])] ^ rpj^g evalua- 
tion of the RHS of Eq. ([M]) can be carried out by the 
Monte Carlo method, similar to that used in ll]. The 
configuration space is a (1 + l)-dimensional N x p lat- 
tice, each direction of which corresponds to space and 
time, respectively, and a spin s* = ±1 is assigned to 
each lattice point. Then the UPO ensemble can be pro- 
duced by the Metropolis algorithm with the "Hamilto- 
nian" pAf[(7A([s]) +/3^([s])], or the content of the brack- 
ets in Eq. (|23p for the Bernoulli CML, by which we can 
obtain the ensemble average and thus A^*. One may 
think that plausible results are not available due to the 
exponentially increasing number of UPOs with N, which 
inevitably restricts the reachable period p to be rather 
short. However, it does not affect since accuracy of the 
Monte Carlo sampling is determined by the proportion of 
dominant orbits in the ensemble average, which depends 
on the product pN (see Eq. (|24p and the form of /([s]) 
below). That is to say, the shortness of sampled orbits 
can be compensated by the large number of DOFs. We 
can therefore compute the GMF by means of the Monte 
Carlo method, over a wide range of {q, (3) through the 
repetition of this step. Note that the step size Ag and 
A/3 must be sufficiently small, otherwise the dominant 
contribution to the average is sustained by the region 
where /([«]) is very small and hence a sampling during a 
practicable Monte Carlo run yields an inadequate result. 
The adequacy of Monte Carlo samplings can be checked 
by plotting H{\ A) by means of Eq. ^} and HH]) and 
seeing that it satisfies the properties of i?(A, A), such as 
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FIG. 4: (Color online) Demonstration of Eqs. (|9]) and (fTO)) for the ID Bernoulli CML with fc = 1, by means of Monte Carlo 
calculations. Lines correspond to (a) mean ^(1,0) and (b) standard deviation —^^(1,0) at A'' = 8,16,32,64 evaluated via 
the GMF ^/(g,/?), which is obtained by averaging results of 400 independent Monte Carlo runs with 100,000 samples after 100 
steps of transients. The range of errors, estimated from standard deviation among the independent runs, is less than 10~^ for (a) 
and 10^"^ for (b), and therefore, negligible. Symbols in both figures indicate results of direct measurement of {j4([s]))upo and 
pNa{A{[s]))'^PQ, respectively, by Monte Carlo simulations with 1, 000, 000 samples after 100 steps of transients. Corresponding 
standard deviations are denoted by error bars. Black dashed curve in (a) represents the exact value (yl)^. Note that plots for 
p = 32 and 64 are nearly at the same place. Lines and symbols for fixed p and are within the range of statistical errors, and 
thus Eqs. (|9]) and (|10[) are confirmed. Furthermore, the two figures show that ^(1,0) and —^^(1,0) converge for sufficiently 
large p, the former of which coincides with the exact value in accordance with Eq. and also for sufficiently large N, which 
indicate that the GMF /3) is analytic in the limit p,N ^ oo. 




FIG. 5: (Color online) The 2nd derivative of the GMF 
with fc = l,p = iV = 16 in the ID BernouUi CML. *(g,/3) 
is obtained from Eq. ^ with Aq = 0.02 and A/3 = 0.01. 
The ensemble average in Eq. ((24} is performed over 50, 000 
Monte Carlo steps after 100 steps discarded as transients. 
The derivative is yielded by the three-point formula. The 
figure shown above is smoothed by taking its moving average 
over 5x5 data points. 



< H{X, A) < X and the concavity. We actually con- 
firmed in Fig. m that the GMF obtained in this way 
satisfies the relations dU and (fTO| . 



C. Analysis of the phase transition 



As is seen from Eq. ^3}, the ID BernouUi CML in- 
volves a 2-dimensional array of spins with the short-range 
interaction in the spatio-temporal configuration space. It 
suggests the occurrence of phase transitions with finite 



values of the parameters. It is indeed the case, which is 
demonstrated by varying q and /3 and plotting the 2nd 
derivative of the "free energy," as is shown in Fig. [5l 
Although we can provide no decisive statement about the 
occurrence of phase transitions from finite-size numerical 
simulations, the two sharp peaks in Fig. [5]clearly indicate 
it, which is confirmed by observing that they grow more 
acute as the system size p and N increases. Therefore 
the ID Bernoulli CML is shown to exhibit phase transi- 
tions in the 2-dimensional space-time. These transitions, 
brought about by varying the temperature parameters in 
the thermodynamic formalism, are called q-phase transi- 
tions in the context of dynamical systems with few DOFs 
[Q]. Moreover the existence of the q-phase transitions 
can be analytically shown in the weak-interaction limit 
fc ^ 0. This can be seen if we neglect O(fc^) in the argu- 
ment of the exponential function in Eq. (|23p. namely 



^ X pNq/2 
V 



n.n. 



(25) 

which results in the canonical partition function of the 2D 
Ising model, where the presence of the 2nd order phase 
transition is certified |29l |. 

Now wc mention the meaning of q-phase transitions ob- 
served in the space-time configuration space in terms of 
the Landau picture of continuous phase transitions. Here 
we do not consider the dependence on a macroscopic ob- 
servable A{x) for the sake of simplicity. The rewriting 
to the thermodynamic formalism with A{x) is straight- 
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forward. First we expand H{X) around the temperature 
q — qa as follows 

H{\) ~Hn + qoX - [B{X - Ao)2 + C(A - Ao)^] , (26) 

where B > and C > because of concavity. The mini- 
mum principle (jl6p yields the Lyapunov exponent at the 
temperature namely X{q) — ^ Xq — (q — qQ)/2B. 
Therefore the 2nd derivative of the "free energy" is 
= 1/2S, which shows that B goes to zero as q ap- 
proaches a 2nd order q-phase transition point. Further, 
the probability to find a positive finite-time Lyapunov 
exponent per 1 DOF A is written as 



(27) 



for large p. Thus the occurrence of q-phase transitions 
involves the breakdown of the central limit theorem for 
finite-time Lyapunov exponents and/or finite-time aver- 
age of macroscopic observable A{x). It can be under- 
stood by the fact that the correlation length and time 
diverge at the 2nd order transition point. Note that, as 
the usual Landau theory, the above statement cannot be 
applied at the vicinity of the transition point 2pNB'^ < C 
[3Q|]. Instead, the occurrence of q-phase transitions is as- 
cribed to the existence of sharp corners in the function 
H{X), which implies an anomaly in the UPO distribu- 
tion with respect to the Lyapunov exponent A ([a;]). Since 
UPOs form the skeleton of the chaotic invariant set 
this means a system accompanying q-phase transitions 
has the invariant set with a global anomalous structure. 
It is worth remarking that q-phase transitions in chaos 
with few DOFs indicate local singularities of the attrac- 
tor, where hyperbolicity is lost 0], whereas q-phase tran- 
sitions in extended chaos treated here signify global ones, 
that are non-analyticities in the distribution function of 
UPOs, which arise without losing hyperbolicity. 

The non-analyticity of the UPO distribution in the ID 
Bernoulli CML can be explicitly confirmed if we consider 
the case k ^ 1 , in which we can refer to the exact solution 
of the 2D Ising model fl^ as is seen in Eq. ([^5]) . In this 
case, a one-to-one correspondence 



Xi[s])^hAi[s]) + Uog2, 



(28) 



reduces H{X,A) to a univariate function H{A). Note 
that a symbol ~ here and in Eqs. (p9)) and (pO)) 
below indicates that both sides of the symbol are 
equal as long as wc neglect O(fc^). Let .^ising 
the partition function of the 2D Ising model per 1 
spin in the thermodynamic limit, namely 2^ising(>^) = 

r 1 

lim„^oo 'I2{sj k}'^^Pi'^J2n.n.'^j,kSj',k') , whcrc n is 
the number of spins. Then we obtain from Eqs. ([7]), p6p . 
(IT7)) . (1^ . and (gS]) the following relations in the limit 
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FIG. 6: The 2nd and 3rd (inset) derivatives of the "entropy 
function" H{A) in the weak-interaction case k <^ 1, which in- 
dicate an anomaly in the UPO distribution of the ID Bernoulli 
CML. The exact solution of the 2D Ising model [29l | is used 
to plot these curves. Note that H{A) does not depend on k 
if we neglect O(fc^). 



p,N ^ oo: 
A{q,f3) 



1 d 



2 dp 

JaW + kq/2), 



logZisi„g(/3-l-fcg/2) 



(29) 



H{A) 0,(^(3+'^^ A +^\og Zi,ing(/3 + kq/2) 

= f],\A)A + i log Zi,i„g(/;, 1(A)), (30) 

where, recalling A([s]) translates into the energy of the 
2D Ising model, the function /a (a;) is monotonic and thus 
its inverse is well defined. The above 2 equations give a 

simple relation = , i-e. the reciprocal spe- 

cific heat. Since the specific heat of the 2D Ising model 
logarithmically diverges at criticality, the 2nd derivative 
of H{A) has a sharp corner and the 3rd derivative di- 
verges, as are shown in Fig. [51 We can of course make the 
same statement on the function H{X) because of Eq. (|28p . 
Therefore the anomalous UPO distribution actually ex- 
ists in the ID Bernoulli CML at least for /c ^ 1, and 
doubtless for all fc, since q-phase transitions are always 
numerically observed. 

On the other hand, as mentioned in Sec. IIII Al the 
CML considered here is equivalent to the ID Ising model, 
so that no Ising phase transition occurs at a finite value of 
the interaction parameter k. It means that, with Eqs. Q 
and Uni), the GMF ^'(g, /3) in the hmit p, N ^ oo is an- 
alytic at the physical situation {q, (3) = (1,0) and the 
CML shows no q-phase transition at that point. Indeed 
Fig. [7] is a phase diagram for several values of k and 
we can see that the transition curves do not go through 
{q,f3) = (1,0) for not so large k. That is, it is true that 
the anomalous part in the UPO distribution exists, but 
at finite k those UPOs are hidden as non-dominant terms 
in the partition function Zi q = X^fxl e^P^^^I^l-* and their 
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FIG. 7; (Color online) q-phase transition curves of the ID Bernoulli CML with p = A'^ = 16 and (a) k = 0.1,0.5, 1.0,2.0, (b) 
k — —0.1, —0.5, —1.0, —2.0, which are indicated by a red solid line, green dot-and-dash line, blue dashed line, and purple dotted 
line, respectively. The black cross is located on the physical situation {q, l3) = (1,0). These transition curves are obtained 
by detecting local maxima of — §^ in q and /3 direction, separately, and then eliminating false maxima which come from 
statistical errors in the Monte Carlo samplings and the finite-size effect. They can be distinguished from rounded-off singularity 
by examining their continuity and dependence on the system size p and N. See also the caption of Fig. [5]for the way to obtain 
The numbers of Monte Carlo samplings are 20, 000, 20, 000, 50, 000, and 100, 000 for \k\ = 0.1, 0.5, 1.0, 2.0 respectively. 



non-analyticity is overwhelmed by the other, analytic and 
dominant terms. However, Fig. [7] shows that, as \k\ is 
increased and goes to infinity, the transition curves move 
and finally reach the physical situation. In other words 
the anomalous UPOs become the dominant terms, and at 
that moment, the Ising phase transition occurs and the 
non-analyticity is uncovered. This is justified by the fact 
that — §^ specifies the lower bound of the fluctuation 
of the Ising energy, or specific heat, so the occurrence of 
q-phase transition at the physical situation just means 
the actual Ising transition. Our consideration reveals the 
role of the anomalous UPO distribution as a "seed" of 
the Ising transition, which is ordinarily hidden. The two 
transition curves, observed at each k in Fig. [71 corre- 
spond to the ferromagnetic (upper curve) and antifer- 
romagnetic (lower curve) transition, respectively, which 
can be understood by comparing transition curves for 
different k in Figs. [7] (a) and (b). 

Finally, we add one comment on the numerical obser- 
vation by Kawasaki and Sasa prj . In order to explain 
the reproduction of macroscopic quantities in turbulence 
from a single UPO 10], they numerically showed that the 
standard deviation of the Ising energy calculated from 
one UPO, i.e. a{A{[x]))upo, goes to zero as the system 
size N increases in the ID Bernoulli CML. This is proved 
by the following facts. Since the model satisfies the scal- 
ing hypothesis (|13p and it does not exhibit a q-phase 
transition at the physical situation for finite values of fc, 
the GMF 5'(q,/3) is assured to be well-defined and ana- 
lytic in the limit p —^ oo and/or N —^ oo. The validity of 
Eqs. ^ and PH)) in both limits is actually suggested by 
means of Monte Carlo calculations as is shown in Fig. |4l 
Therefore, by taking the limit N ^ oo in Eq. (jlOp . we 
obtain a(A([a;]))upo even for a finite period p. This 
is what Kawasaki and Sasa numerically observed 



and might be a ground for the macroscopic reproduc- 
tion in turbulence [To| . In other words, any hyperbolic 
extended systems which satisfy Eq. (|13p , or the large de- 
viation theorem, possess this property. We can also see 
from Eq. (jlOp and Fig. HJ^b) that the accuracy of a sin- 
gle UPO estimate, i.e. standard deviation a{A{[x]))upo, 
asymptotically scales as {pN)~^^^. Note that, however, 
the period must not be too short (p > 32 in the case 
of the ID BernouUi CML) in order to regard the UPO 
average (A([a;]))upo as a good approximation of the tur- 
bulent average (A)^ (see Fig. |3]and Eq. 



IV. ANALYSIS OF THE ID REPELLING CML 
-SOLVABLE CASE- 

The UPO expansion and the thermodynamic formal- 
ism dealt with in Sec. [Til are also applicable to repelling 
systems insofar as we concentrate our attention into the 
dynamics on chaotic invariant sets. A modification is 
required only on Eq. which is replaced by 

lim *(1,0) = a, (31) 

p — ^oo 

where a is the escape rate per 1 DOE of the repeller, 
i.e. Zi,o = Ewe"^^^^[^l^ e-^^". The space-time 
"Hamiltonian" pNX([x]) can be now constructed at will 
without the strong constraint of Eq. ^ , hence a solvable 
model is available. 

Here we adopt a ID coupled repeller map lattice in- 
troduced by Just and Schmiiser [17|. Dynamical vari- 
ables are x* £ I = [—1,1] defined at each site i = 
0,1,-- - ,N — 1, with a "spin" variable 

s* = sgnx*. (32) 
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FIG. 8: The local map, given by Eqs. (|33j and ((Ml, for the 
repelling CML by Just and Schmiiser. ^t^^ indicates a 

slope of each piecewise linear part. 



The time evolution of x* is yielded by a piecewise linear 
map 



(33) 



with 



/ \ Jr — .s{x + c) if-l<a;<-c, 
I 1 -+.s{x + c) it — c < x < 0, 



, , _ I s{x — c) if < X < c, 

i ++.s{x - c) it c < a; < 1, 



(34b) 



for some values of slopes r±±_s > 1, as functions of 
the neighboring spin s, and a constant < c < 1. The 
periodic boundary condition = Sg is considered to 
close the definition. The local map defined in this way is 
sketched in Fig. [S] The invariant set of this CML can be 
symbolized again via the partition 



Vs ^ g-'s, (I) X 9-]s, (/) X • • • X g;^_^^,^ (/) , (35) 



where the local partition gj^^i+ii-^) depicted in Fig. 
IS] Hence the first and second symbol in the subscripts 
of r±±_s indicate a spin of site i at time t and t + 1 
respectively. 

Now we apply the thermodynamic formalism to the 
model. The slopes F^t^t+i , which prescribe both the 
local dynamics and the interaction between neighboring 
sites, can be chosen arbitrarily provided that the local 
map does not cross the boundary x = 0, ±1. Here we 
choose the simplest form after Just and Schmiiser (l7j . 



exp [—J si (s 



t+i 



eo 



(36) 



where J and eo are some constants. A macroscopic quan- 
tity is set to be the Ising energy again, namely Eq. P^ . 
Thus, the partition function ^ for this model is calcu- 
lated as 



{W-lp-l ^ 
1=0 t=0 J 



(37) 



r 



which is nothing but the canonical partition function for 
the 2D Ising model on the square lattice with anisotropic 
interaction. Note that the Bernoulli CML treated in the 
previous section results in the 2D Ising model only at the 
weak interaction limit, whereas for the repelling CML it 
holds for all q and (3. By setting (q,/3) = (1,0), i.e. 
physical situation, in Eq. ([57)) . the model turns out to 
show the 2D Ising transition at J = Jc = ^ In (l-f V2) as 
is shown in [13]. Moreover, the existence and location of 
the q-phase transition curve (ffe, 0c) is also given exactly 
by Onsager's celebrated paper [2^ as 

/3c = ^ Arcsinh [1 / sinh(2(?c J)] - QcJ- (38) 
Figure [5] shows its phase diagram for several values of the 



coupling constant J. In this case, the transition curve 
passes the physical situation (g, /3) = (1,0) linearly as 
J goes through Jc, which can be explicitly written from 
Eq. §^ as 

for fixed Qc, (39a) 
for fixed Pc- (39b) 



/3e^-2(J- J,), 

gc =i 1 - (J - Jc)/Jc 



This linear dependence of the q-phase transition point 
{qc, f3c) on the control parameter J at the vicinity of the 
actual phase transition point J = Jc results from the fact 
that the transition observed here is not a marginal one. 
Hence this relation between the two transition points is 
expected to be general. 

On the anomalous structure of the invariant set with 
respect to the UPO distribution and its role in the oc- 
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q 

FIG. 9: (Color online) q-phase transition curves (|38p of the 
ID repelling CML defined by Eqs. ((32l)-(l34l) and JMl), with 
J = 0.3, Jc ~ 0.44, 0.6, 1.0, which are drawn with a red solid 
line, green dot-and-dash line, blue dashed line, and purple 
dotted line, respectively. The black cross indicates the physi- 
cal situation {q, (3) = (1, 0). The transition curves for negative 
J are obtained by reflecting the figure over q- or /3-axis. 

currence of the phase transition, the same statement as 
the previous section holds, which can be demonstrated 
directly for this model since the rigorous solution is avail- 
able. 



V. DISCUSSION 

The intimate relation between q-phase transitions and 
phase transitions in the sense of statistical mechanics is 
investigated on the basis of the thermodynamic formal- 
ism and the UPO expansion. Since mathematically the 
partition function ^ has the identical form to that of 
the canonical statistical mechanics with q and /3 as in- 
verse temperature, many useful relations in equilibrium 
physics, such as Eqs. ([9]), (fTO| . and (fT7|) . also hold in 
extended chaotic systems, which can be far from equi- 
librium. Although similar relations have been already 
pointed out for dynamical systems with few DOFs by 
several authors 0, ^ ^ , we reconstructed it for extended 
systems with concerning the number of DOFs N explic- 
itly. By that means the analogy is kept with the equi- 
librium statistical mechanics of several-dimensional sys- 
tems. A richer harvest may be reaped from it, as long 
as attention is paid to the strict constraint of Eq. ([8]) for 
systems without escape. 

As regards phase transitions, anomalously distributed 
UPOs turn out to be responsible, which show sharp cor- 
ners in the distribution fimction and can be visualized in 
terms of q-phase transition. The anomalous part exists 
in systems with transitions, over the range of control pa- 
rameters where the topological structure of an attractor 
does not change. It is ordinarily hidden as non-dominant 
terms in the partition function Zi^ and no critical behav- 
ior is observed there. The actual transition occurs when 
the control parameters are varied, the UPO distribution 



is changed and finally the anomalous UPOs become a 
dominant part. In this sense we call the anomalous part 
of the UPO distribution "seed" of phase transitions. 

One question may arise here: "What brings this 
anomalous UPO distribution to dynamical systems with 
phase transitions?" The answer is clear for the two Ising- 
like systems considered in this paper, where the origin of 
phase transitions is by construction well known from the 
knowledge of equilibrium statistical mechanics: the com- 
petition between interaction energy and entropy is rel- 
evant. Taking ferromagnet for example, the free energy 
F{E) — E — TS{E) calculated under some fixed energy is 
increased by low entropy for strongly ferromagnetic con- 
figurations (corresponding to low E) , while it is raised by 
high internal energy for strongly paramagnetic configu- 
rations (corresponding to high E). It means there are 
intermediate configurations where the two mechanisms 
compete. In fact, this competition occurs at one point, 
i.e. at some specific value of E, in the thermodynamic 
limit, which brings a sharp corner to the functional form 
of S{E). Then the phase transition occurs at a temper- 
ature which minimizes the free energy at that point. 

The role of UPOs in g-phase transitions — not in ac- 
tual phase transitions — is exactly the same as that 
of microstates which we have seen above. That is, the 
anomalous distribution of UPOs results from the compe- 
tition between average positive Lyapunov exponent and 
topological entropy. This mechanism may be widespread 
even among "natural" extended chaotic systems, because 
it is reasonable to expect that the number of UPOs with 
plenty of large Lyapunov exponents is very small, and 
that it grows in the same manner as equilibrium mi- 
crostates (see Eq. ([13]) and Ref. [l[3l|). Note that the 
existence of symbolic dynamics is also not required, since 
the underlying basis described in Sec. |lT] is constructed 
generally for hyperbolic maps. 

Taking into account the above considerations and the 
aforesaid similarities in statistics of macroscopic observ- 
ables such as Eqs. ^ and (fTU)) . we can propose a def- 
inition of phase transitions and their order in extended 
dynamical systems in the Ehrenfest's sense: phase tran- 
sitions are associated with the singularity of the GMF 
at the physical situation {q,f3) = (1,0). The transition 
can be said to be of n-th order if an n-th derivative of 
the GMF with respect to g or /3 does not exist or has 
a discontinuity, and hence the system is accompanied 
by n-th order q-phase transition. This is a mathemat- 
ically simple-minded statement as well as that based on 
the non-uniqueness of a natural measure. Moreover it is 
worth remarking that the proposed definition can clearly 
characterize both first and higher order transitions, while 
definitions in the Gibbsian sense have some ambiguity 
when it comes to treat higher order transitions. Though 
the new definition has also several problems which will 
be mentioned later, it can be used to classify phase tran- 
sitions out of equilibrium and to investigate their nature 
further. This is the main proposition of this paper. 

The second outcome is on the propriety of treating lo- 
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cally defined control parameters as macroscopic "temper- 
ature." As we have seen in Sec. IIIII and lTVl systems which 
exhibit a phase transition have an anomalous UPO dis- 
tribution, the position of whose non-analytic part is spec- 
ified by critical nominal temperatures qc and f3c. There- 
fore, our observation that they vary with local control pa- 
rameters means that a change in local parameters leads 
to a change in "macroscopic temperature" through the 
non-analyticity of the UPO distribution. This macro- 
scopic temperature actually takes part in phase tran- 
sitions, since it crosses through the physical situation 
= (IjO) when a transition occurs, and since it 
mathematically works in the same way as real tempera- 
ture in equilibrium systems (cf. Eqs. ([6]), (l9|), and (fTO|) ). 

Let us then discuss the replacement of "temperature" 
by local control parameters around transition points. Let 
P denote some control parameter in an extended dynam- 
ical system with a phase transition. As we have already 
seen in Sec. IIVI the relation between parameter P and q- 
phase transition point {qc, (3c) at the vicinity of the tran- 
sition point P = Pc is expected to be linearly dependent 

P-f3c^CiiP- Pc), for fixed qc, (40a) 
q^qc^C2{P-Pc), for fixed Pc , (40b) 

for transitions which are not marginal. Here we set q 
and P at the physical situation {q, /3) = (1, 0). Therefore, 
as far as some universal relation in equilibrium physics is 
concerned which is not affected by microscopic details of 
models, e.g. critical behavior, the same relation may hold 
in extended chaotic systems by replacing inverse temper- 
ature {l/ksT - l/ksTc) with Ci{P - Pc)- A similar 
statement could be also said on the positive Lyapunov 
exponent per 1 DOF, in which case Eq. I|40b[) is used for 
the replacement. Conditions which should be satisfied at 
least by the relation are that (1) it is about some macro- 
scopic quantities obtained by differentiating a free energy, 
and that (2) its mathematical expression itself is insensi- 
tive to variations in the control parameter P. The differ- 
ence in the relation between 2nd or higher order moments 
and derivatives of the "free energy" from its counterpart 
in equilibrium statistical mechanics might also have an 
influence — in extended chaos the derivatives can only 
tell the lower bound of the corresponding moments due to 
temporal correlation, as is seen in Eq. pI7|) — . Provided 
that those restrictions are taken into consideration, we 
believe that the mentioned replacement can be applied 
to a wide range of extended systems. Note that univer- 
sal scaling relations in critical behavior satisfy the above 
conditions and thus corresponding critical exponents are 
likely to be kept invariant under the replacement of tem- 
perature T by the local control parameter P. It can be a 
basis on which scaling relations indeed work under such 
a replacement in some high-dimensional chaotic systems 
(e.g. [l,[il). 

In order to justify the above arguments on a rigid basis, 
several problems need to be clarified. To begin with, it 
is unclear in extended chaotic systems how common the 



existence of the anomalous UPO distribution is, and also 
how prevalent its relation to phase transitions is. The 
latter is especially important, since it involves a change 
in the UPO distribution and thus there is no counterpart 
in equilibrium statistical mechanics. Further studies are 
crucial. 

From a fundamental point of view, we do not math- 
ematically care in this paper either the existence of the 
two limits p — > cxD and N ~f oo, their order or the fact 
that they do not commute. They are undoubtedly im- 
portant in order to argue spatio-teniporal chaos on the 
mathematically proper basis 0, IH, [H, [H, [2l| . An ex- 
amination of the behavior of infinite-size systems requires 
that we first take the limit iV — » oo and then p — > cx), 
at variance with usual statistical mechanics where the 
limit is taken over sizes of all dimensions simultaneously. 
This may be the reason why some extended chaotic sys- 
tems defined in c?-dimensional space show critical behav- 
iors of d-dimensional universality classes despite the cor- 
responding {d + l)-dimensional configuration space p^ . 
The problem of the incommutability should be consid- 
ered seriously, since the definition of phase transitions by 
means of the singularity of the GMF involves both the 
limit p, N —^ oo. 

Another problem is on the arbitrariness for the choice 
of a macroscopic quantity A{x) when we deal with q- 
phase transitions with respect to /?. There is no clear 
criterion for it, except for that A{x) must be affected 
by phase transitions : its expectation value or fluctua- 
tion must show a discontinuity or divergence. An order 
parameter of the considered transition is a candidate. In 
our examples in Sec. Illll and HVl however, we adopted the 
quantity which can be regarded as energy on the anal- 
ogy with equilibrium spin systems. The relation between 
the choice for A{x) and the behavior of q-phase transi- 
tion curves near critical points remains to be clarifled, 
especially about their linear dependence on control pa- 
rameters such as Eq. (PH)) . 

In conclusion, the old concept of the thermodynamic 
formalism and the periodic orbit expansion turns out to 
be useful to characterize phase transitions in extended 
dynamical systems. Theoretically, one possible definition 
of phase transitions is proposed, which is complemen- 
tary to the usual definition in terms of a change in the 
number of natural measures. It can be used to classify 
and examine non-equilibrium transitions in chaotic sys- 
tems, especially higher order transitions, with the help 
of a suitable technique to generate or approximate the 
UPO ensemble. Recently developed method such as Ref. 
[3^.[33l| might be applied for this purpose. With regard to 
experiments and numerical simulations, a ground is ob- 
tained on which we can sometimes treat an externally im- 
posed control parameter as macroscopic "temperature" 
around phase transition points. Although some signifi- 
cant problems are left for future studies, this assertion is 
expected to support discussions on universality classes in 
non-equilibrium systems from real and numerical exper- 
iments. 



13 



Acknowledgments with K. Kaneko, M. Kawasaki, and S. Tatsunii. One 

of us (K. T.) would also like to thank K. Nakajima for 
The authors gratefully acknowledge sincere comments letting him use the PC cluster Cenju for this work, 
from H. Tasaki and S. Sasa, and beneficial discussions 



[1] D. Ruelle, Thermodynamic formalism (Addison- Wesley, 
1978). 

[2] C. Beck and F. Schlogl, Thermodynamics of chaotic 
systems: an introduction (Cambridge University Press, 
1993). 

[3] T. C. Halsey, M. H. Jensen, L. P. Kadanoff, I. Procaccia, 
and B. I. Shraiman, Phys. Rev. A 33, 1141 (1986). 

[4] H. Fujisaka and M. Inoue, Prog. Tiicor. Piiys. 77, 1334 
(1987); Piiys. Rev. A 39, 1376 (1989); Phys. Rev. A 41, 
5302 (1990). 

[5] P. Cvitanovic, Phys. Rev. Lett. 61, 2729 (1988); 

P.Cvitanovic et al, Chaos : Classical and Quantum 

(Niels Bohr Institute, Copenhagen, 2005), available at 

http : // chaosbook . org/. 
[6] H. Hata, T. Horita, H. Mori, T. Morita, and K. Tomita, 

Prog. Theor. Phys. 80, 809 (1988); Prog. Thoer. Phys. 

81, 11 (1989); T. Horita, H. Hata, H. Mori, T. Morita, K. 

Tomita, S. Kuroki, and H. Okamoto, Prog. Theor. Phys. 

80, 793 (1988); S. Sato and K. Honda, Phys. Rev. A 42, 

3233 (1990). 

[7] See e.g. L. A. Bunimovich and Ya. G. Sinai, Nonlinearity 
1, 491 (1988); for hyperbolic maps, e.g. M. Jiang, Nonlin- 
earity 8, 631 (1995); for reviews, e.g. L. A. Bunimovich, 
Physica D 103, 1 (1997); J. Bricmont and A. Kupiainen, 
Physica D 103, 18 (1997). 

[8] See, e.g. P. Manncvillc, in Macroscopic Modelling of Tur- 
bulent Flows, Lecture Notes in Physics Vol. 230, edited by 
U. Frisch, J. B. Keller, G. Papanicolaou, and O. Piron- 
neau (Springer- Verlag, Berlin, 1985), p. 319; F. Chris- 
tiansen, P. Cvitanovic, and V. Putkaradze, Nonlinearity 
10, 55 (1997). 

[9] A. Politi and A. Torcini, Phys. Rev. Lett. 69, 3421 
(1992). 

[10] G. Kawahara and S. Kida, J. Fluid Mcch. 449, 291 
(2001); S. Kato and M. Yamada, Phys. Rev. E 68, 
025302(R) (2003); L. van Veen, S. Kida, and G. Kawa- 
hara, Fluid Dyn. Res. 38, 19 (2006). 

[11] M. Kawasaki and S. Sasa, Phys. Rev. E 72, 037202 
(2005). 

[12] H. Chate and P. Manneville, Europliys. Lett., 17, 291 
(1992); Prog. Theor. Phys. 87, 1 (1992); H. Chate, A. 
Lcmaitrc, P. Marcq, and P. Manneville, Physica A 224, 

447 (1996). 

[13] P. Marcq, H. Chate, and P. Manneville, Prog. Theor. 

Phys. SuppL, 161, 244 (2006). 
[14] J. Miller and D. A. Huse, Phys. Rev. E 48, 2528 (1993). 
[15] P. Marcq, H. Chate, and P. Manneville, Phys. Rev. Lett. 

77, 4003 (1996); Phys. Rev. E 55, 2606 (1997). 
[16] G. Gielis and R. S. MacKay, Nonlinearity 13, 867 (2000); 



R. S. MacKay, in Dynamics of coupled map lattices 
and related spatially extended systems, Lecture Notes in 
Physics Vol. 671, edited by J. -R. Chazottcs and B. Fer- 
nandez (Springer- Verlag, Berlin, 2005), p. 65. 

[17] W. Just, J. Stat. Phys. 105, 133 (2001); W. Just and F. 
Schmiiser, in Dynamics of coupled map lattices and re- 
lated spatially extended systems, Lecture Notes in Physics 
Vol. 671, edited by J. -R. Chazottes and B. Fernandez 
(Springer- Verlag, Berlin, 2005), p. 33. 

[18] M. Blank and L. A. Bunimovich, Nonlinearity 16, 387 

(2003) . 

[19] J. -B. Bardet and G. Keller, Nonhnearity 19, 2193 
(2006). 

[20] A. C. D. van Enter, R. Fernandez, and A. D. Sokal, J. 
Stat. Phys. 72, 879 (1993). See Sec. 2.6.5 and references 
therein. 

[21] J. Bricmont and A. Kupiainen, Commun. Math. Phys. 

178, 703 (1996). 
[22] J. L. Lebowitz, C. Maes, and E. R. Speer, J. Stat. Phys. 

59, 117 (1990). 

[23] P. C. Hohcnbcrg and B. I. Shraiman, Physica D 37, 109 
(1989); M. S. Bourzutschky and M. C. Cross, Chaos 2, 
173 (1992); F. Sastre, I. Dornic, and H. Chate, Phys. 
Rev. Lett. 91, 267205 (2003). 

[24] T.Kai and K. Tomita, Prog. Theor. Phys. 64, 1532 
(1980); T. Morita, H. Hata, H. Mori, T. Horita, and K. 
Tomita, Prog. Theor. Phys. 79, 296 (1988); C. Grcbogi, 
E. Ott, and J. A. Yorke, Phys. Rev. A 37, 1711 (1988). 

[25] M. Kawasaki (private communication). 

[26] M. Sano, S. Sato, and Y. Sawada, Prog. Theor. Phys. 76, 
945 (1986); J. -P. Eckmann and I. Procaccia, Phys. Rev. 
A 34, 659 (1986). 

[27] H. Sakaguchi, Prog. Theor. Phys. 80, 7 (1988). 

[28] D. Frenkel, in Molecular-Dynamics Simulation of 
Statistical- Mechanical Systems, Proceedings of the Inter- 
national School of Physics "Enrico Fermi", Course 97, 
edited by G. Ciccotti and W. G. Hoover (North-Holland, 
Amsterdam, 1986), p. 151. 

[29] L. Onsagcr, Phys. Rev. 65, 117 (1944). 

[30] L. D. Landau and E. M. Lifshitz, Statistical Physics, 3rd 
ed. (Pergamon Press, Oxford, 1980). 

[31] R. Kubo, Statistical Mechanics: an advanced course 
with problems and solutions (North-Holland, Amster- 
dam, 1965). 

[32] Y. Lan and P. Cvitanovic, Phys. Rev. E 69, 016217 

(2004) . 

[33] S. Sasa and K. Hayashi, Europhys. Lett. 74, 156 (2006). 



